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Presentation  Outline 


A)  Review  of  Soil  Model  governing  equations 

B)  Development  of  pedo-transfer  functions  (terrain 
database  to  engineering  properties) 

C)  Ethan  Allen  Firing  Range  (EAFR)  terrain  database 
conversion 

D)  Dynamic  Soil  Model  example  (Stryker  tire  with  EAFR 
soil) 
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A)  Soil  Model  Governing  Equations 


1)  Boussinesq/Cerruti  vertical  stress  distribution 

2)  Soil  response  -  Vertical  visco-elasto-plastic 
displacement 

3)  Shear  stress/deformation 

Slip  (longitudinal),  Turning  (lateral) 

Bulldozing 
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Disturbed  width  for 
tracked  vehicles 
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Compression  Model  Approaches 


"Essentially,  all  models  are  wrong,  but 
some  are  useful/'  —  George  Box 


Wheel  numerics  -  Dimensionless  numbers 
Bekker  equation 

Boussinesq  stress  and  Compression  Index 
FEM  or  DEM 
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Boussinesq  Vertical  Stress  Distribution 
due  to  Vertical  Force 


3  Pz3 

2  x[(x2  +  y1  +  z2)f 2 
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from  Feda,  1978 


Boussinesq  Stress  Distribution 


•  Assumes  soil  follows  the  theory  of  elasticity 

•  Assumes  contact  area  of  the  tire  is  a  rectangular  plate 

•  Modified  equation  with  Frohlich's  concentration  factor 


vWzv 

2;t(j c2+j2+z2)(v/2+1) 


(Source:  Feda,  1978) 


LOAD :  7.5  KNl  1690  Ub 


(Source:  Wong,  2008) 


2)  Soil  response  -  Vertical  visco-elasto- 

plastic  displacement 

•  Compression  index 

•  Influence  of  soil  type  and  moisture 

•  Rebound 

•  Time  dependent  sinkage 
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Clog {vj°k) 


re  factors 


10 

-2 


from  Larson  et  al.,  1980 


Elasto-plastic  deformation  (rebound) 
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from  Upadyhaya  et  al.,  2002 


3)  Shear  stress/deformation 

Lateral/longitudinal  stress-strain,  shear 
deformation  modulus  (Janosi/Wong) 

Bulldozing  (soil  displacement)  analysis,  passive 
lateral  earth  pressure 
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Lateral  shear  deformation  modulus 

slip  displacement 

=  rmax  •  (l  -  e~J/K  )  =  (c  +  a  •  tan  <f) ■  (l  -  e~J',K  ) 

•  Where  ris  the  shear  stress, 

•  j  is  the  shear  displacement, 

•  c  is  the  internal  cohesion  of  the  soil, 

•  ^  is  the  angle  of  internal  friction  of  the  soil, 

•  K  is  defined  as  the  shear  deformation  modulus 
(Wong,  2001). 
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"Bulldozer"  forces 

(passive  lateral  earth  pressure) 


P  =  b('A  yZ2N ,  +  2  czJn^) 

N^=  tan2  (45  +  ^/ 2) 
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B)  Development  of  pedo-transfer  functions 


Engineering  parameters  needed  by  soil  model 

-  compression  index 

-  rebound  index 

-  time  constant 

-  cohesion 

-  friction  angle 

-  shear  deformation  modulus 

-  density 
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Ethan  Allen  Firing  Range  Database 
(what  we  are  starting  with) 


Elevation  >=  0 


Hex  code 

Decimal 

code 

Terrain 

Condition 

Soil  Type 

Moisture 

content 

code 

RCI  0-6  psi  RCI6-12 

Snow 

Depth 

Snow 

Density 

Frost 

Depth 

Thaw 

Depth 

4A00 

18944 

NFG 

ML 

NOR 

50 

80 

0 

0 

0 

0 

4A02 

18946 

NFG 

ML 

NOR 

50 

80 

0 

0 

0 

0 

Elevation  >  152 

Hex  code 

Decimal 

code 

Terrain 

Condition 

Soil  Type 

Moisture 

content 

code 

RCI  0-6  psi  RCI  6-12 

Snow 

Depth 

Snow 

Density 

Frost 

Depth 

Thaw 

Depth 

5300 

44308 

NFG 

CL 

SLP 

80 

80 

0 

0 

0 

0 

5300 

21248 

NFG 

CL 

SLP 

80 

80 

0 

0 

0 

0 

5304 

21252 

NFG 

CL 

SLP 

80 

80 

0 

0 

0 

0 

530C 

21260 

NFG 

CL 

SLP 

80 

80 

0 

0 

0 

0 

5314 

21268 

NFG 

CL 

SLP 

80 

80 

0 

0 

0 

0 

532C 

21292 

NFG 

CL 

SLP 

80 

80 

0 

0 

0 

0 

5334 

21300 

NFG 

CL 

SLP 

80 

80 

0 

0 

0 

0 

5340 

21312 

NFG 

CL 

SLP 

80 

80 

0 

0 

0 

0 

5344 

21316 

NFG 

CL 

SLP 

80 

80 

0 

0 

0 

0 

Terrain  database  descriptors 
(off-road  soil  conditions): 

Terrain  Condition  -  Fine-grained  soil  NFG,  coarse¬ 
grained  soil  NCG,  road,  other. 

Soil  Type  -  Uses  Unified  Soil  Classification  System 
(USCS)  CL,  CH,  ML) 

Moisture  content  code  -  NOR  (normal)  and  SLP 
(slippery)  for  fine-grained  soil,  normal  for  coarse¬ 
grained  soil) 

RCI  0-6  -  Soil  strength  0-6  inches,  RCI  for  fine¬ 
grained  soils,  Cl  for  coarse-grained  soils. 

RCI  6-12  -  Soil  strength  6-12  inches,  RCI  for  fine¬ 
grained  soils,  Cl  for  coarse-grained  soils. 
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Pedo-transfer  function 


Need  to  transfer  existing  terrain  database 
descriptors  to  soil  engineering  properties. 
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Unified  Soil  Classification  System  (USCS) 


MAJOR  DIVISIONS 

GROUP 

SYMBOLS 

TYPICAL  NAMES 

FIELD  INDENTIFICATION  PROCEDURES 
(excluding  particles  larger  than  3  inches  and 
basing  fractions  on  estimated  weights) 

INFORMATION  REQUIRED 
FOR  DESCRIBING  SOILS 

1 

2 

3 

4 

5 

6 

Coarse-grained  Soils 

More  than  half  of  material  Is  larger  than 

No  200  sieve  size 
matlest  particle  visible  to  the  naked  eye. 

♦  £ 
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ll  tl 
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is 

(Clean 
Gravels 
Little  or 
no  fines) 

GW 

Weil-graded  gravels,  gravel-sand 
mixtures,  little  or  no  fines 

Wide  range  in  grain  sizes  and  substantial  amounts  of  all 
intermediate  particle  sizes 

For  undisturbed  soils  add  information 
on  stratification,  degree  of 
compactness,  cementation,  moisture 
conditions,  and  drainage 
characteristics. 

Give  typical  name:  Indicate 
approximate  percentage  of  sand  and 
gravel  maximum  size,  angutanty, 
surface  condition  and  hardness  of 
the  coarse  grains;  local  or  geologic 
name  and  other  pertinent  descriptive 
information,  and  symbol  in 
parentheses 

Example: 

Siltv  sand  aravellv:  about  20%  hard, 
angular  gravel  particles  1/2in. 
maximum  size;  rounded  and 
subangular  sand  grains,  coarse  to 
fine;  about  15%  non  plastic  fines  with 
low  dry  strength;  well  compacted  and 
moist  in  place;  alluvial  sand  (SM). 

GP 

Poorly  graded  gravels  or  gravel-sand 
mixtures,  little  or  no  fines 

Predommady  one  size  or  a  range  of  sizes  with  some 
intermediate  sizes  missing 

(Gravels  with 
Fines 

Appreciable 
amounl  of 
fines) 

GM 

SUty  gra  vels,  gravei-sand-silt  mixtures 

Nonplastic  fines  or  fines  with  low  plasticity  (for 
identification  procedures  see  ML  below) 

GC 

Clayey  gravels,  gravel-sand-clay  mixtures 

Plastic  fines  (for  identification  see  CL  be! cm/) 

Sands 

More  than  half  of 
coarse  fraction  is 
smaller  than  No.  4  sieve 
size 

(For  visual  classification, 
eqiiv8lent  to 

Clean 
Sands 
(Little 
or  no 
Ines) 

SW 

Well-graded  sands,  gravelly  sands,  little  or 
no  fines 

Wide  range  in  grain  sizes  and  substantial  amounts  of  all 
intermediate  sizes  missing 

SP 

Poorly  graded  sands  or  gravelly  sands, 
little  or  no  fines 

Predominantly  one  size  or  a  range  of  sizes  with  some 
intermediate  sizes  missing 

Sands  with 
Fines 

(Appreciable 
amounl  of 
fines) 

SM 

Silty  sands,  sand-silt  mixtures 

Nonplastic  fines  or  fines  with  low  plasticity  (for 
idenbficabon 

SC 

Clayey  sands,  sand-clay  mixtures 

Plastic  fines  (for  idenbficabon  procedures  see  CL  below) 

Fine-grained  Soite 

More  than  half  of  material  is  smaller  than 

No  200  sieve  size. 

I  he  No.  200  sieve  size  Is  about  the  s 

Identification  Procedures 
on  Fraction  smaller  than  No.  40  Sieve  Size 

Dry  Strength 
(Crushing 
Characteristics) 

Dilatancy 
(Reaction  to 
Shakina) 

Toughness 
(Consistency 
near  PL} 

Silts  and  Clays 

Liquid  limit  is  less 
than  50. 

ML 

Inorganic  silts  and  very  fine  sands,  rack 
flour  silty  or  dayey  fine  sands  or  clayey 
silts  with  sliqht  plasticity 

None  to  slight 

Quick  to  slow 

None 

For  undisturbed  soils  add  information 
on  structure,  stratification 
consistency  in  undisturbed  and 
remolded  states,  moisture  and 
drainage  conditions. 

Give  typical  name,  indicate  degree 
and  character  or  plasticity,  amount 
and  maximum  size  of  coarse  grains, 
color  in  wet  conditions,  odor  (if  any), 
local  or  geologic  name,  and  other 
pertinent  descriptive  information,  and 
symbol  in  parentheses. 

Example: 

Clavev  silt,  brown,  sllahtiv  olastic: 
small  percentage  of  fine  sand; 
numerous  vertical  root  holes;  firm 
and  dry  and  place  loess  (ML). 

CL 

Inorganic  clays  of  low  to  medium  plasticity 
gravelly  clays,  sandy  days,  silty  days, 
lean  days 

Medium  to  higfi 

Norte  to  very  slow 

Medium 

OL 

Organic  silts  and  organic  silty  days  of  low 
plash  city 

Sl  ight  to  medium 

Stow 

Slight 

Soils  and  Clays 

Liquid  limit  is 
greater  than  50. 

MH 

Inorganic  silts,  micaceous  or 
diatomaceous  fine  sandy  or  silty  soils, 
elastic  silts 

Slight  to  medium 

Slow  to  none 

Slight  to 
medium 

CH 

Inorganic  days  of  high  plasticity,  fat  days 

High  to  very  high 

Nate 

High 

OH 

Organic  days  and  silts  of  medium  to  high 
plasticity 

Medium  to  high 

None  to  very  stow 

Slight  to 
medium 

Highly  Organic  Soils 

Pt 

Peat  and  other  highly  °rg^r)i|5J^>jl^\SS  I F 1 E 

Readily  identified  by  color,  odor,  spongy  feel  and 

Lrequentty  by  fibrous  texture 
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The  centroid  percent  Sand,  Silt,  and  Clay  and 
particle  size  for  9  USCS  Classes 


USCS  Soil 

Type 

■ 

%Sand 

Avg 

Average  Grain 

Size  (mm) 

Clay  Index  (0-5) 

Grain  Size 
Index  (0-5) 

CH 

73.3 

13.3 

13.3 

0.35 

n 

CL 

47.2 

24.0 

28.9 

0.76 

3.21 

2.32 

MH 

34.8 

50.0 

15.2 

0.41 

2.38 

4.27 

ML 

13.6 

61.5 

24.9 

0.67 

0.93 

2.63 

SC 

26.2 

8.8 

65.0 

1.69 

1.79 

1.04 

SM 

8.8 

26.2 

65.0 

1.70 

0.60 

1.04 

SW/SP-SC 

6.9 

2.3 

90.8 

2.36 

0.47 

0.75 

SW/SP-SM 

2.3 

6.9 

90.8 

2.36 

0.16 

0.75 

SW&SP 

1.7 

1.7 

96.7 

2.51 

0.11 

0.70 
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Ethan  Allen  Firing  Range  terrain  database  over  1000  soil  units  (rows). 
Ignoring  the  elevation  and  non-soil  units,  a  set  of  13  unique  soil  units 
were  identified  as  shown  below. 


Pedo-transfer  functions 


(p  =  24.03-3.06  •  {Clay  Index)  + 1.03  •  {Grain  Size  Index)  +  .0047  •  {RCI) 
c  -  8.65  + 12.15  •  {Clay  _  Index)  -  8.94  •  {Grain  Size  Index)  +  .0070  •  {RCI) 

C  =  .2175  +  .01245  •  {%Clay)  -  .000131  •  ( %Clay  f 

K  =  2 1 .74  -  3 .22  •  {Clay  _  Index)  +  0.47  •  ( GSI ) 

=  1.68-0.10- {Clay  _  Index )  -  0.04  •  {GSI) 

ReboundConstant  =  0.0 11  +  0.004  •  (GSI) 

pK  =  1.544  -  0.00556  •  {%Clay)  -  3.468  x  10"5  •  {%Clayf 

ST  =0.003461  + 1.742  xlO"4  -(%Silt)-  2.980x1  O'6  -{%Silt)2  (C<M) 

ST  =  0.003217  +  3.251xl0“4 -{%Clay)~  5.385  xlO"6 -{%Clay)2  (s) 
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C)  Predicted  Soil  Engineering  Properties  for  the  EAFR  Terrain  Codes 


Terrain 

Soil 

Type 

Moisture 

RCI 

Clay 

Index 

Grain 

Size 

Index 

Friction 
Angle  (2) 

NFG 

CL 

SLP 

80 

3.21 

2.32 

17.0 

NFG 

GM 

NOR 

150 

l 

l 

NFG 

ML 

NOR 

100 

0.93 

2.63 

24.4 

NFG 

ML 

NOR 

65 

2.63 

24.2 

NFG 

ML 

NOR 

43.5 

m 

2.63 

24.1 

NFG 

ML 

NOR 

30 

m 

2.63 

24.0 

OTH 

PT 

NOR 

17.5 

■ 

■ 

1 

1 

1 

1 

NFG 

SM 

NOR 

130 

1.04 

23.9 

NFG 

SM 

NOR 

110 

m 

1.04 

23.8 

NFG 

SM 

NOR 

65 

0.60 

1.04 

23.6 

NCG 

SP 

NOR 

130 

■H 

m 

NFG 

CL 

SLP 

130 

3.21 

2.32 

17.2 

NFG 

SC 

SLP 

80 

1.79 

1.04 

4.6 


27.8 


21.6 


Bulk  Rebound 

Density,p  Constant 

C(g/cm3)  (g/cm3) _ (g/cm3) 
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D)  Dynamic  Soil  Compression  Model 


•  A 1  cm3  soil  element,  5  cm  below  the  soil 
surface,  positioned  along  the  tire  centerline 

•  EAFR  terrain  -  CL  (USCS),  80  psi  (RCI), 
SLP  (moisture) 

•  Simple  Stryker  tire  (22.25  kN  normal  load) 
traveling  at  3  m/s  with  a  tire  inflation 
pressure  of  276  kPa  (40  psi).  (uniformly 
distributed  rectangular  surface  load) 


EAFR  Soil  Engineering  Properties 


•  Terrain  -  NFG 

•  Soil  Type  -  CL 

•  Moisture  -  SLP 

•  RCI  -  80 


•  Clay  Index  -  3.21 

•  Grain  Size  Index  -  2.32 

•  Friction  Angle  (9)  - 17.0 

•  c  (kPa)  -  27.4 

•  K  (cm)  - 10 

•  C  (g/cm3)  -  0.51 

•  Bulk  Density,  p  (g/cm3)  - 1.25 

•  Rebound  Constant  (g/cm3)  -  0.020 

•  Time  Constant  (s)  -  0.2 
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350 
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Stryker  Soil  Stress  at  0.05  m  -  Centerline  -  (3  m/s) 
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Soil  Element  Vertical  Displacement  (cm) 
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Bulk  Density  (g/crrt3) 
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Model  E  Output 
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Model  P  Output 
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Soil  Element  Vertical  Displacement  (cm) 


Number  of  Wheel  Passes 


0 1.5  m/s  Travel  Speed  03.0  m/s  Travel  Speed  06.0  m/s  Travel  Speed 
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Change  in  Energy  (J) 


Number  of  Wheel  Passes 


01.5  m/s  T ravel  Speed  03.0  m/s  T ravel  Speed  06.0  m/s  T ravel  Speed 
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Peak  Power  (W) 


Number  of  Wheel  Passes 


El  1 .5  m/s  Travel  Speed  El 3.0  m/s  Travel  Speed  D6.0  m/s  Travel  Speed 
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Travel  Speed  (m/s) 


BTotal  Vertical  Displacement  (cm)  B Total  Energy  (J)  BPeak  Power  (W) 
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Benefits  of  New  Soil  Model 


•Theoretical  (soil  deformation,  energy,  power) 

•  Not  dependent  of  scale 

•  Real-time  modeling 

•  Include  visco-elastic  (time  constant/rebound) 

•  Uses  existing  terrain  database  with  pedo- 
transfer  functions. 


Questions 


**Disclaimer:  Reference  herein  to  any  specific  commercial  company,  product,  process  or  service  by 
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endorsement,  recommendation,  or  favoring  by  the  United  States  Government  or  the  Department  of 
the  Army  (DoA).  The  opinions  of  the  authors  expressed  herein  do  not  necessarily  state  or  reflect  those 
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